Accepted for publication in AJ 

Preprint typeset using I^'T^]X style emulateapj v. 04/20/08 



LOCATING THE YOUNGEST H ii REGIONS IN M82 WITH 7 MM CONTINUUM MAPS 
Chao-Wei Tsai\ Jean L. Turner\ Sara C. BECK^ David S. Meier^-"', & Paul T. P. Ho''-'* 

Draft version March 10, 2009 

ABSTRACT 

We present 7 mm Very Large Array continuum images of the starburst galaxy M82. On arcsecond 
scales, two-thirds of the 7mm continuum consists of free-free emission from Hii regions. In the 
subarcsecond resolution map, we identify 14 compact sources, including 9 bright H ii regions with 
Nhyc > 10^^ sec~^ . Four of the Hii regions have rising spectra, implying emission measures > 
10® cm~^ pc. Except for one compact source with peculiar features, all other compact radio sources 
are found in dust lanes and do not have optical or near-infrared continuum counterparts. Four regions 
of extended, high brightness (EM > 10^ cm~^pc) radio emission are found in our high resolution 
map, including some as large as ~2", or 30 pc, representing either associations of small Hii regions, 
or sheetlike structures of denser gas. The good correlation between 7 mm emission and Spitzer IRAC 
8 fim continuum-removed PAH feature suggests that PAH emission may track the recently formed 
OB stars. We find an excellent correlation between molecular gas and star formation, particularly 
dense gas traced by HCN, down to the ^ 45 pc scale in M82. We also find star formation efficiencies 
(SFEs) of 1-10% on the same scale, based on CO maps. The highest SFE are found in regions with 
the highest dense gas fractions. 

Subject headings: Hli regions - galaxies: individual (M 82, NGC 3034) - galaxies: ISM - galaxies: 
starburst - galaxies: star clusters - radio continuum: galaxies 



1. INTRODUCTION 

Free-free emission from H ii regions at radio wave- 
lengths is a good tracer of young and heavily embed- 
ded regions of massive star formation such as luminous 
starbursts and nuclear star forming region s, where vi- 
sual extinctions can be tens ( Ho et al.|1990 1 to hundreds 
of magnitudes, affecting even near and mid-infrared ob- 
servations. Extragalactic continuum emission is compli- 
cated by the presence of non-thermal synchrotron emis- 
sion within the beam, part icularly at wavelengths longer 



than 6 cm (Condon 19921. While it is possible in the- 
ory to do a spectral and spatial separation of free- free 
and synchr otron emission ( Turner fc Ho|1983 1994 Mc- 
Donald et al. 2002) with niulti- wavelength observations, 
the free-free emission can be most easily measured at 
millimeter wavelengths, where radio synchrotron has de- 
clined and dust has not yet begun to take over. The 
enhanced high frequency capabilities of the Very Large 
Array (VLA) — improved K and Q band receivers, fast- 
switching calibration techniques — allow us to detect 
and resolve the free- free emission from large H ii regions 
in nearby galaxies. 

We report VLA observations of the 7 mm continuum 
emission of the starburst galaxy M82. An interacting 
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dwarf with a bar and gas ring in the central kpc ( Telesco 
et al.|1 991 Lar kin et al.|1994| [Ixhtermann fc Lacy|1995| 
Shen fc Lo,1995 1, M82 i s a pro digious former of stars, at a 
rate of ~ 10 Mq yr'^ ( [Teles co 1988). The burst of star 
formation may have been triggered by the interaction 
with M81 ( jYun et al.ll993|, or by infalling m olecular gas 
driven by the stellar oar iG 



]Greve et al. 20021. Its current 



active star form ation can be traced back to 50 Myr ago 
(de Grijs 20011. The strong star formation activity is 
probably responsible for the hot gas outflows fo und in 
optical, mm, and X-ray studies (O'ConncU & Mangano 
19781 INakai et al. 1987: Loiseau et al. 1990: Matsumoto 



et al 



20001 IMatsushita et al. 



2000j .Stevens et al.„2003| ). 



Our goal is to use the 7 mm maps to isolate compact 
Hii regions in M82, and to determine their free-free flux 
density, from which we can infer Lyman continuum rates, 
or N^yc- With the subarcsecond resolution afforded by 
the VLA we can determine locations and sizes of the 
bright Hii regions on sca les of a few pc. E xtinction is 



high and patchy in M82 
to be Ay 



Simon et al 



1979k, estimated 



50 (Forster ISchreiber et al. 2001 1 on large 



scales, and as muc h as several hund red mag behind the 



molecular clouds ( Keto et al.||2005 l, thus radio free- free 
emission is a particulaxly valuable tracer of ionized gas 
in this dusty starburst galaxy. 

Spectral energy distributions (SEDs) of the compact 
radio sources at longer, cm wavelengths (Wills et al. 
1997 1 lAllen fc Kronberg] [T9981 [Menl|l999] |McDonald 



et a. 



l2C)02) in M82 have shown them to be mostly 
non-thermal synchrotron sources, probably supernova 
remnants (SNR ) but possibly synchrotron-em itting wind 
driven bubbles ( |Seaquist fc Stankovic||2007 ). The struc- 
tures and the expansion properties ot these non-thermal 
sources have been revealed by Very Long Baseline Inter- 
ferometry (VLBI) and Multi-Element Radio Linked In- 
terferometry (MERLIN) with 3-35 mas resolutions and 



2 



Tsai et al. 



multi-epoch monitoring at 18 and 6 cm ([P edlar et 



McDonald et al. 2001 Beswick et al 



al 

I'tenech 



2008 1. In addition to the radio continuum work 



1999] 
et al. 

Rodrigucz -'Rico et al.| ( r2004[ ) studied the H92q; and H52q! 
radio recombination hues m M82, including the 7 mm 
(42.952 GHz) continuum, in 0.6" or ~ 10 pc resolution. 

In this paper, we have 7 mm images at two resolutions: 
high resolution, 0.2", for compact structures on scales of 
3 pc, and low resolution, 1.5", for structure on scales of 
> 25 pc. The map with 0.2" synthesized beam is the 
highest resolution map of M82 at millimeter wavelengths 
so far. Lower resolution maps are important for judging 
the amount of missing, undersampled emission in higher 
resolution VLA images, which act as spatial high-pass 
filters. 

The distance to M82 was recently determined to be 
3.6±0.3 Mpc from the optical study of Cepheids in M81 
using HST (Frecdman et al. 1994), assuming the same 
distance to both galaxies. Observations of the tip of the 
red giant branch in M82 itself suggests 3.9 ±0.3 (sys) 
±0.3 (random) Mpc ( [Sakai fc Madore||1999| . For con- 
sistency with previous work, we adopt the 3.6 Mpc dis- 
tance, with 1"=17 pc. At this distance a 1 mJy thermal 
source at 7 mm represents an H ii region ionized by the 
equivalent of ^ 150 07 stars. 

2. OBSERVATIONS 

The radio data were acquired with the NRAO VLA^ 
using A- and B- configurations at 43 GHz (NRAO 
program ID: AT303) and D-configuration at 45 GHz 
(AM839, PI: Meier, D. S.). Weather during the A-array 
observations in October 2004 was not good enough for 
useful 7 mm work. 

The "high resolution" (sub-arcsecond) maps we 
present in this paper were made in the B-configuration 
on 22 April 2005 in continuum observation mode (effec- 
tive bandwidth ~172 MHz). The calibration cycle time 
for fast-switching was set to 100-120 seconds. 3C147 
was used for calibration of the absolute flux scale and 
0958+655 was the phase calibrator. Uncertainty in the 
absolute flux density scale is estimated to be <5%, based 
on VLA estimates of the bootstrap accuracy at 7 mm. 
For most of our sources, statistical uncertainty due to 
noise dominates the uncertainty in flux density scale. To- 
tal on source time is 1.3 hours. The {u,v) sampling covers 
30-1610 fcA, resulting in a ~0.2" beam, the size of which 
varies slightly with weighting. For our B-array map, the 
estimated 0maxi the largest angular scale of structures 
that our observation is sensitive to, is about 3.5". 

Lower resolution 7 mm continuum observations (hence- 
forth "low resolution maps") were made with the VLA 
D-array at 45.47 GHz on November 27 2005 in spectral 
line mode, as part of a spectral line program (AM839). 
The setup was two 25 MHz spectrometer tunings per IF 
set side-by-side. The uv dataset was generated by fitting 
the uv spectrometer data minus the edge channel with a 
first order polynomial using UVLSF task. The net band- 
width is ^ 31 MHz. Because IF 2 has contamination 
from a recombination line, only IF 1 data is included. 
The observations used fast- switching with 180 seconds 
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cycle time. 3C147 and 3C286 were used as flux density 
calibrators, and 0958±655 was the phase calibrator. To- 
tal integration time on source is 2.2 hours. The sampled 
{u,v) plane covers 5-157 kX, resulting in a ~1.5" beam. 
We estimate the O^ax is ~ 20" in our D-array map. 

The maps were calibrated and reduced in AiPS follow- 
ing high frequency reduction procedures, and were de- 
convolved with CLEAN. The noise level in the maps is 
~ 0.1 mJy bm~^ for the high resolution (B-configuration) 
maps and ~ 0.5 mJybm"^ for the lower resolution (D- 
configuration) maps, but varies slightly depending on 
weighting. Flux measurements and positions were ob- 
tained with IMFIT task in AiPS unless otherwise noted. 
The AIPS task PBCOR was used to correct the flux den- 
sity scale of these data for the primary beam efficiency 
of the instrument. 

3. OVERVIEW: COMPACT AND DIFFUSE EMISSION 
SOURCES AT 7 MM IN M82 

The high-resolution B-array map (Fig. Ilk) covers 
the inner 700 pc x 350 pc (40" x 20") of M82 with a 
0.19" X 0.15" (~ 3 pc) beam. Structures larger than 
3.5" are suppressed in this map. We identify 14 compact 
sources with flux densities > 0.5 mJy (Tft=12 K; 5 a), 
distributed in a ~ 40" x 10" east-west s wath with posi- 



tion a ngle along the optical major axis (Kronberg et al. 



1972|. The compact sources are high brightness, dense, 
"Tu 



and luminous H ii regions and young supernova remnants 
(SNR). For the Hii regions, our 5cr detection limit corre- 
sponds to N^yc ^ 0.6 X lO'^^ (assuming optically thin 
free-free emission), corresponding to about 80 O stars, or 
a cluster such as R1 36. The compact o ptical H ii region 
M82-A1 studied by [Smith et al] ( [2006[ ) falls just below 
our detection limit. 

Properties of the compact sources in natural weighted 
map are g iven in Table [ l l We adopt the source designa- 
tions of , Kronberg et arj^l985 ), using the B1950 coordi- 
nate of each source relative to 9''51'"+69°54', to facilitate 
comparisons with earlier studies, but the J2000 coordi- 
nates are given also. In Table |2] we have collected the 
radio continuum flux densities of all 14 compact sources 
from the literature (18 - 1.3 cm Unger et al. 



nberg et al.||1985 
McDonald et al 



Allen fc KronbergI 



2002P and this work. The flux den 



]fT984llKro] 
Allen 199^ 



sity measurements at different wavelengths are done at 
similar resolutions. The beam sizes of maps at 18 cm to 
1.3 cm range from 0.2" to 0.3". The 7 mm measurements 
quoted in Table |2] (as well as data shown in Figure J2] [3| 
are made on a tapered map of 0.2" resolution in order to 
match beam sizes at other wavelengths from the litera- 
ture. The quoted numbers in Table [2] are total integrated 
flux densities of compact sources, and the effect of dif- 
ferent beam is less than 30%. In Figure [2] [3] and Table 
[2] we show the SEDs of all 14 compact sources at wave- 
length 18 cm to 7 mm. The details of model fitting are 
included in Table [3) Our discussion of properties of these 
compact sources, including their SEDs, can be found in 
pi] and fp] 

The high-resolution B-array map also reveals four 
regions of comparatively extended emission (EM = 
J n^dl). The bright, high frequency free-free emission 
detected here arises from gas with very high emission 
measure, so it is surprising to find such extended sources. 
They are 1-2" across, or 17-35 pc. They are similar in 
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extent to the large H ii region associated with NGC 3603 
in our Galaxy, but with a much higher brightness and 
EM. We discuss the signi ficance of these high brightness, 
extended regions in §4.3| 

The low resolution D-array map is shown in Figure 
This map, with a beam of 1'.'99 xl'.'47 p.a. 17°, 
is sensitive to lower brightness emission, Tf, ~ 0.1 K 
rms (cr '^O.S mJy/bm) and also more extended emission 
structures (up to 60", covering the entire field shown) 
than the high resolution map. The low resolution im- 
age shows four major peaks with intensity 12 mJy/bm 
or higher. The integrated flux density in the naturally 
weighted D-array map is 1.2 ± 0.1 Jy. The total flux re- 
ported for M 82 from WMAP at 41 GHz flux density is 
1.3 ± 0.2 Jy ( [Bennett et"al]|2003| . Thus our low resolu- 
tion D array map recovers, within the uncertainties, all 
the 7 mm flux density. The diffuse emission structure is 
elongated in the same direction and over the same ex- 
tent as the swath of compact sources, and agrees with 
the orientation of the stellar ell iptical component s een in 
the deep near infrared images ( [Mayya et al][2005| . 

We have fltted 17 peaks in the low resolution map 
with Gaussians, and the corresponding flux densities, 
locations, and the sizes of the peaks are listed in Ta- 
ble |4 All of these pea ks are found in the longer wave- 



The 14 compact sources of the high resolution map 
(Fig.jlk) contribute < 10% of the total 7 mm flux density 
in thelow resolution D-array map (Fig.jlJ)). These com- 
pact sources in the 7 mm high resolution B-array map 
also appear in lower frequency radio cont inuum maps 
dMuxlow et~ar]|1994| [H uang et aLlll994| [Wills ct aL] 
l007| [Allen fc Kronberg 1998 Allen 1999: McDon^Tet] 



20bl[ [Fenech et al.[2008J. We have computed^ radio 



Shvbs for the 14 sources, wnich are shown in Figures 
and[3| with details of the SED modeling given in Tab 
3| We identify nine regions with spectral index a > —0.2 
(X v"') as Hii regions, and the others as SNR. The 
flux densities at 43 GHz and the spectral indices are 
consis tent with previous SEDs and source character iza- 



and 
Ol re- 



length VLA images of Alleli[ ( 1999 1 except peaks W4a 
and El a, which are newly identified in this work. The 
peaks sum to ~ 0.65 Jy at 7 mm. If this flux density 
associated with low resolution peaks is all optically thin 
free-free emission, and Tg — 10,000 K, the inferred Ly- 
man continuum photon rate (NLyc) is ~ 8.4 x 10^^ sec~~^. 
(For Te of 5,000, NLyc would increase by 30 ~ 40%.) 
This value is in good agreement with that estimated 
from t he fre e-free continuum at longer wavelengths by 



tions (Alle n fc Kronberg|1998l [McDonald et al.[2002] ) ex- 
cept for the two non-thermal sources 43.99 2-1-59. ?[) anc 
41.942-1-57.62, which are discussed in §6.1[ and 
spectively. 

Of the nine H ii regions in the high resolution map, two 
have flux densities at more than five frequencies. Fit- 
ting of the spectra of these two sources suggests that 
they have free-free turnover frequencies (where r ~ 1) 
Vt ^^-"2.5-5 GHz. To have f( ~ a few GHz requires emis- 
sion measures of 2 — 9 x 10^ pc cm~^. 

The other sources have flux densities at only three fre- 
quencies, and a single power law represents their ob- 
served flux densities fairly well (Fig. |2|. Some of the 
deviations may be due to small amounts of extended 
emission in the sources and the slightly (15-20%) differ- 
ing beams of the flux densities from the literature. Four 
sources have rising SED (a > 0). Rising spectrum radio 
sources have been found in many extragalactic systems, 
and are thought to be giant compact Hii regions with 



high electron density (Turner et al. 1998). These are 



wavelength spectral de composition modeling), a nd at 



shor ter wavelengths b y Car lstrom & Kronberg ( 1991 1 
and [Matsushita et al.j ( [2005| ), who find 6.9 x lO'^^ sec^ 
and 7.9 x 10'^'^ sec~^ at 3 mm and 2.6 mm, respec- 
tively (these numbers are interferometric and therefore 
are lower limits). It also agr ees with that sugge sted 
from recombination line studies ( Seaquist et al. 1985 ) and 
the extinction-c o rrecte d mid-infrared fiuxes of 11* orsterl 
(20011. 



Schreiber et al. 



Allen (1999 9.1 x 10 sec^ , from 1.3 - 20 cm multi- much larger than Galactic compact H ii regions because 

the super star clusters (SSCs), the clusters exciting them, 
are much larger both in excitation and in physical size. 
For the sources here, the emission appears to remain par- 
tially optically thick at frequencies up to 43 GHz, so we 
expect to see brightness temperatures of 10^ — 10'^ K. 
Their observed brightness temperatures are lower than 
this, and imply EMs > 3.5 x 10^ pc cm~^ (derived from 
peak fiux density), consistent with turnover frequencies 
of 3 GHz, an order of magnitude lower than observed. 
The low brightness temperatures indicate small filling 
factors, and signal the presence of very compact sources 
of significantly less than 3 pc extent, which we cannot 
resolve. 

Considering the radio continuum luminosity, bright- 
ness of the emission, and the turnover frequency, these 
appear to be young H ii regions, and we infer that the star 
clusters exciting these nebulae are also relatively young. 
The total flux density of just the detected compact H ii 
regions comes to ~ 38 mJy, which is ~ 6% of the total 
thermal flux density at 7 mm (at the same sensitivity) , in 
other words, ~ 6% of the ionizing flux density. But the 
high resolution 7 mm image fails to detect all of the com- 
pact Hii regions in M 82, even fa irly optically luminous 
ones such as M82-A1 ([Smith & Gallagher 2001), so we 



The good agreement in N^yc 
between our analysis and previous studies indicates that 
the peaks we see in the low resolution map of Figure [1]d 
dominate the total free-free emission. 

Non-thermal emission at 43 GHz is weak, but con- 
tributes a non-negligible fraction of the integrated flux 
density in M82. The estimate of the non-thermal compo- 
nent is 35% of the total at 7 mm, extrapolating from VLA 
cm- wave flux densities of the same region (Allen 1999) . 
This would give a thermal flux density of ~ 0.8 Jy, 15 
- 20% higher than 0.65 Jy we obtain from summing the 
peaks in the D-array map of Figure Wp. We attribute 
the discrepancy to thermal emission from some H ii re- 
gions distributed on too large a spatial scale to be seen 
in peaks in these images. 



STAR FORMATION TRACER: 7 MM CONTINUUM 
EMISSION IN M82 

4.1. The Compact Thermal Radio Sources 



believe that there are more compact regions below our 
current detection limit. In this case the Hii region phase 
in M82 must last at least ~ 6% of the ~10 Myr lifetime 
of O stars, and possibly longer, to a Myr or more. To 
restrain the Hii regions from expanding requires exter- 
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nal or internal confinements, perhaps external high pres- 
sures (de Pree et al. 1995 Smith et al. 2006), or in the 



larger cl usters, potentially slowing by gravity ( Turner et 
aL pOOSl ) 



4.2. Lack of Optical/Near-IR Continuum Counterparts 

The compact sources are shown in Figure |4] overlaid 
on a color (VRI) image from the Spitzer Legacy project 
SINGS (Kennicutt ct al. 2003 ). Nearly all of the compact 
sources are tound m dustTKhes. The only exception is 
source 43.922+59.70, which is located near the dynamical 
center, and a potential AGN candidate (§6.) 

We have attempted to align the compact radio sources 
identified in this and previous work (Allen 1999, McDon- 
ald et al 2002, and references therein) wit h the published 
HST n ear-I R continuum images of |Alonso-Herrero et al.| 
( |2003| and |McCrady~et al. (2003|, the catalogued op- 
tical clusters ( [Melo et al.||2005| ), and the archived HST 
optical measurements, 'i'here is no consistent alignment 
of the radio and other sources within the HST pointing 
uncertainty. The H ii regions we see are presumably as- 
sociated with young SSCs; that the clusters are not seen 
in the optical/Near-IR continuum argues that the extinc- 
tion must be high {Ak > 6) . The molecular clouds alone 
can contribute Ay of 10^-10^ if they happen to be along 
the line of sight to the H ii regions. 

If the average extinction at arcsecond resolu tion is 
Ay 14 (from H-K color by |McLeod et al.|1993 assum- 
ing Av/Ak = 10), the line-ot-sight extinctions implied 
in the lack of near-IR sources to these young clusters 
must be 4 times higher. This is consistent with the high 
average {Ak ^ 1) and patchy extinction derived on large 
(11", or ~ 200 pc) scales toward th e H ii regions i n M82 



Rieke et al. 



(20031, 



from Brack ett line observations by Willner et al. ( 1977 
Ay = 25), [Simon et alTI (|1979| Ay = 14) 

( |1980| Ay = 26), and [Forster Schreiber 

Ay ~ 50), much of which appears to be internal to the 
H II regions themselves. 

In the Pa a image of Alonso-Herrero et al. 
we identify 3 (39.227+54.32, 40.933+58.9S, and 
42.072+58.52) out of our 9 compact thermal sources, 
and one non-thermal source (43.992+59.70) as associ- 
ated with compact Pa a structures. If the Pa a emission 
is associated with the corresponding radio Hii regions, 
we estimate, based on the inferred Lyman continuum flux 
densities, that the extinction toward these IR- visible neb- 
ulae is Ak ^ 9. None of these three sources are apparent 
in the HST near-IR continuum images, and it agrees with 
the high extinction suggested in the Pa a image. 

4.3. High EM Extended Structures at 7 mm: Sheets of 

Hot Gas 

Besides the compact dense sources in Table |3] we de- 
tect extended regions in M82 with relatively high EM. 
In Figure [5] we show B-array maps of free-free emission 
regions with integrated flux density of 50 - 70 mJy in 
5"- 9" area. These regions have EM = 1.2 - 2.0 x 
lO'^ pc cm~^ . This agrees with the D-array measure- 
ments shown in Table U 

High EM can be produced either by high rms elec- 
tron density or by a long sightline {dl) of ionized gas. If 
high electron density is responsible, then to reconcile the 
high EM with the observed flux {N^yc) requires that the 



line of sight length be much smaller than the observed 
transverse distance. For example, to maintain the ioniza- 
tion of a spherical volume of radius 15 pc, the observed 
transverse diameter of source in Fig |5|i, and rms elec- 
tron density 10^ cm~^ (> 5 x 10'^ cm was inferred by 



Rodriguez- Rico et al.|[2004, in a 2" resolution) requires a 



Lyman continuum rate of N^yc ^ 2 x 10 s . The ob- 
served flux of this structure indicates iV^yc ~ 10^^ s~^. 
This implies that the actual emitting volume is smaller, 
i.e., the emitting region is shcetlike. There are a number 
of such patches evident in Figure [Sj While these large 
regions of high electron density are dynamically stable 
in the high pressure P/k > 10^ cm~'^ K environment of 
the M82 starburst, the excitation requirements are in- 
consistent with other star formation indicators, with the 
possible exception of the "yin-yang" shaped source (Fig. 
[5fc) located at the peak of the low resolution map (Fig.lTl. 
.hese structures could be the ionized surfaces of chim- 



neys ([Wills et al. 1999b I or superbubbles ([Neininger et 
lal.'l9^8"Wei]j ei al.lii>yi?[ [Wills et al.|1999bl[MaisusTu!a 

jCt al. 2000 2005 1 . The coincidence of their locations with 
dust lane and GO gas distribution implies that they are 
associated with in the molecular clouds, which would be 
consistent with an origin in the outer, ionized layers of 
large molecular clouds. 

Despite its location near the edge of the superbubble, 
the "yin-yang" source is an exception in that it does ap- 
pear to have sufficient flux to support the io nization, 
and there is a bright mid-IR source here to o (Achter- 
mann fc Lacy[[l995 Lipscy & Plavchan 2004). However 



this source has an unusual structure that is coherent and 
bubble-like on a scale of 1.5" = 26 pc. If this were a bub- 
ble, an energy of 10^°-^ ergs would be required to clear 
the inner region of gas. 

Another explanation for these extended regions is that 
they represent a collection of smaller Hii regions rep- 
resenting distributed star formation at the end of the 
stellar bar in M82. Regions (b) - (d) in Fi gure [S] are 
distribut ed on the 10" ionized ring of M82 (Larkm et 
al. 19941 [Achtermann fc Lacylll995| , which 



, , , , n iay repre- 
sent the a:2-orbits ot a bar ([i talcsco ct al. 11991 Larkin et 
al [ [1994j [Achtermann fc Lacy||1995, jGreve et al.||2002p 



Tne region (a) is on the east ionized ridge, and probably 
near the end of the X2-orbits. Star formation is observed 
to be near bar ends, perhaps facilitated by cloud-cloud 
collisions. 

4.4. Hn Regions, Molecular Cas, and Dust 

4.4.1. Correlation with Molecular Gas: CO(l-O) and 
HON (1-0) Line Emission 

The dense star forming regions traced by high bright- 
ness 7 mm free-free emission are likely to be still close to 
their natal giant molecular clouds (GMGs) . We therefore 
expect a high correlation of thermal free-free emission 
and molecular line emission, although at some sufficiently 
high resolution, they should begin to anti-correlate. In 
Figure |6] we show the 7 mm map with contours of C0(1- 
0) line emission from the OVRO map of Shen fc Lo ( 1995 1 
and an HGN(l-O) map from the Plateau de Bure Interter- 
ometer (Schilke, private communication). The molecular 
gas tracer GO(l-O) is more extended in the plane of the 
galaxy, while the dense molecular gas tracer HGN(l-O) 
emission is confined to the star formation area. It is 
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apparent that the correlation of the 7 mm continuum 
emission, largely free-free emission, with HCN gas is ex- 
cellent. 

A quantitative comparison of the molecular line fluxes 
and the free free flux conflrms this conclusion. Figure 
[T] shows the beam-to-beam comparison between the low 
resolution 7 mm emission and HCN. The linear correla- 
tion between HCN and radio continuum shown in Fig- 
ure [7] agrees with the tight correlation between global 
HCN emission and star formation tracers (such as far-IR 
emission ) di scovered in local normal spiral and starburst 
galaxies [Cao fc Sol omon (20 04a| , in Galactic molecular 
dense cores J Wu efal.,2005J , and in high redshift galaxies 
( Gao et al.||2UU7p . Uur plot shows that this good corre- 
lation between dense molecular gas and star formation 
extends down to scales of tens of pc. 

To m ore directly compare our results with [Gao et aT] 
(20071) who determined the correlation in terms of in- 



frared luminosity, L/jj and HCN line flux, we have trans- 
lated the 7 mm flux density into L/^. We assume an 
0-star-dominated stellar population which yields the re- 
lation of LiR ^ Lqb ~ 2 X 105 Lq X (iTpfe^), fol- 
lowing [Turner et aT] (|2000| and [Beck et al.] p002[ ), or 
LiR ~ 2.9 X 10^ Lq X (^^) at 3.6 Mpc. This relation 
could result in underestimates if evolved stellar popula- 
tions contribute signiflcantl y to Lir. Our 7 mm values 



are plotted with the data of Gao et al. ( 2007 1 on Figure 
[8j Further discussion of this correlation is m |5.2[ 



4.4.2. Correlation with Spitzer 8 /im Emission 

The polycyclic aromatic hydrocarbon (PAH) emission 
features in the mid-infrared are expected to trace photo- 
dissociation regions (PDRs) gene rated by the softer UV 



photo ns produced by OB stars (Hollenbach & Tielens 
|1997) . Emission in the 8 /im band of the IRA C cameraTs 
oo 



ominated by aromatic feature emission (Engelbracht 
al. 20061, probably PAHs. The contribution of starlig 
to the 8 fim band is estimat ed to be minimal , at ^10% 
(based on IRAC 3.6 /Ltm flux; [Helou et al |2004[ ) . The res- 
olution of IRAC at this wavelength, 1.7 , is comparable 
to our low resolution 7 mm image. 

M82 was observed at 8 fim by IRAC on Spitzer as part 
of the Spitzer Infrared Nearby Galaxy Survey (SINGS; 



Kennicutt et al. 2003 ) . The 7 mm contours of our low 
resolution U-array map are shown overlaid on the IRAC 
8 /im map in Figure [9] The 7 mm and 8 /im emission 
have very similar morphologies. 

We have compared the corrected flux density in IRAC 
8/im and at VLA D-array 7 mm beam-by-beam. The 
IRAC 8/im flux density is corrected to remo ve the stellar 
emissi on in the IRAC 8/zm band, following Helou et al.| 
(|2004|. The 7 mm map is convolved to 2.2" x 2.2" to 



match the PSF of IRAC channel 4 {8^J.m). The hfim v.s. 
Ijmm plot is shown in Figure [Tol We limit the data points 
shown to those with 3 a detection at 7 mm. This process 
reduces the contribution of non-thermal radio emission 
which is fainter at high frequency and more extended. 
The ratio is linear, with Isnm/ Iimm ~ 45 and a Isfim 
zero intercept ~ 90 ±10 mJy/bm in 2.2" beam. The high 
correlation of 8 /im and 7 mm emission may mean that 
they are excited by the same relatively strong radiation 
sources on the ~ 40 pc length scale. The /s/im that 
appears even at zero Irmm, about 40% of the total, could 



be excited by the relatively older population of B stars 
that do not produce strong radio continuum emission. 

Integrating the correlated (above the offset level) 8 /im 
flux density gives about 30 Jy. Our beam-to-beam linear 
correlation between /g/xm and Irmm suggests that the 
Lyman continuum rate is 



Lyc 



2.2 X 10^^ 



sec 



'Mpc 



correlated 



-1 ^ ^ ^2^'^8/im 



Jy 



for Isfj,m > 0.06 mJy/pc^. Other relations between 
Ssiim and Niyc have been found from work at different 
wavelengths. Deriving a Paa flux from our Niyc (using 
Niyc = 9.8 X 10^3 gg^- 



eras sec~ 



from 



Ho et al. 



1990 and Storey fc Hummer||1995| ), the flux de nsity ex- 
pected froni our 7 mm derived I\lvc using the Calzetti 
et al. (20071 relation for Ssfim/ Spaa would give about 
a tactor ot o more S's/zm than observed. This may not 
be out side the large scatter in the observed relationship, 
which Calzetti et al. (20071 discuss. Another prediction 
of the 8/im flux density can be made from the empiri- 
cal relation s N^yd Lq,2PAh ~ 10^® sec~^ /Lq (Peeters 



|et al.|[2004[ based on studies on normal galaxies, star 
burst galaxies, and Hii r egions in the Loca l Group) and 
^6.2 nm/Ls.6 /im = 1-5 ( [Smith et alT||2007[ ). This would 
give S's/im about a factor ot iU lower than observed. All 
relations of the PAH luminosity to the Niyc are uncertain 
and have large scatter because the PAH flux will depend 
on the geometry of the radiation field, the role of dust, 
the presence of B stars, and possibly other factors not 
yet considered. The spatial correlation of the PAH fea- 
ture and the ionization in M82 is clear; the quantitative 
relation is not yet so. 

PAH emission may be as sociated with a B star popu- 
lation (Peeters et al. 2004) because the strong UV field 
of O stars could destroy PAHs. If we constrain the ion- 
ization sources of the radio continuum emission in H ii 
regions to be B stars, forcing the initial mass function 
(IMF) to be deficient in O stars, then to produce the ob- 
served Nj^yc would require the total mass of young stars 
to increase by 2 orders of magnitude relative to a normal 
IMF. This leads to a star formation rate of few hundred 
solar mass per year for the past 100 Myr, which is highly 
unlikely. 

Why do the radio continuum emission peaks correlate 
with 8 fim peaks, where the PAH carriers should be de- 
stroyed by the UV? It is possible that dust within the 
H II regions around the massiv e stars absorbs a signif- 
icant fraction of the UV flux (Forster Schreiber et al. 
20011. Or the production rate of PAHs may be higher 
than the destruction. To examine the later scenario, we 
estimate the lower limit to the number density of hydro- 
gen atoms, uh, for PAH survival in the UV photon fields 
observed: 



Uh > FFUV,5.f)-l3.6eV /{H X c) 



(2) 



where c is the speed of light, Fpuv is the photon fiux rate 
[photon cm~^/sec] for photons with energy between 5- 
13.6 eV, and the dimensionless H is defined as the Habing 
photodissociation parameter*, following §12.7 of Dopita 



° represents the ratio of local FpuY photo density and the 
density of atomic hydrogen. 
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fc Sutherland! ( |2003[ ). [Dopita fc Sutherland! ( [2003[ ) sug- 
gest H ~ (J. (Jo, an upper limit trom l ine-ot-sight average 



from observation fAUai n et al.||1996 1 , for PAH molecule 
destruction. The Starburst99 code (with Salpeter IMF) 
gives a ratio of the far UV ionizing flux, Npjjy , to Niy^ 
of Npijv /Niyc ~ 3 for 0-type dominated stellar popula- 
tion, thus we find 



Fpijy[cm ^ sec 



1.3 X 10' 



mjy/hm 



(3) 



So the lower limit of hydrogen atom number density, 
uh = 900 X ( "^T^/T ) cm~^ for PAH survival in the lo- 
calized region where It mm is observed. In M82, the peak 
7 mm VLA flux density in a ~ 2" beam ~ 5 — 20 rnJy, 
requiring uh > 5 x 10^ cm~^ . li H ^ 0.005 which en- 
sures the PAH survival in the PDRs but not inside the 
ionized media (Do pita et al.^ 2005), nn > 5 x 10** cm~^ 
is required to supply the observed PAH, close to the crit- 
ical density for HCN excitation. This agrees with the 
good correlation between intensities of PAH with HCN 
emission in M82 found in 2.4" beam (not shown in this 
paper). 

5. STAR FORMATION HISTORY AND EFFICIENCY 
5.1. Star Formation and its Recent History in M82 



Forster Schreiber et al. ( 2003 1 proposed the star forma- 
tion history ot selected st ar burst regions in M82 based on 
ISO mid-IR imaging and spectroscopy for A = 1-45 /xm. 
For the central 16" x 10" region, they found star forma- 
tion had peaked near 5 and 9 Myr ago with amplitudes 
of 6.3 and 18.5 Mq yr-' 
to less than 0.1 yr 



^, respectively, and had dipped 
^ at 6 Myr or so. Previously 



Forster Schreiber et al. (20011 had commented on the 
lack of the Wolt-Kayet ( VV-K) sjDCctral signatures in M82, 
which may be due to high extinction toward the "young" 
star forming regions (3-6 Myr) where the W-R stars live, 
but could also reflect the star formation rate (SFR) dip 
near 3-1- Myr ago. Our observations suggest that high 
extinction is the reason for the absence of WR features, 
since all of the young radio free-free sources are found in 
heavily visually obscured regions. 
The integrated flux density at 7 mm of the full region 



Forster Schreiber et al. ( |2001[|2003| mapped is about 340 
mjy in the low resolution 7 mm map. This equates to a 
total mass in young stars of 1.3 x lO'' Mq formed within 
the average lifetime of O stars of 1-3 Myr, for a Salpeter 
IMF of 1-100 Mq. The contribution of the previous two 
starburst events at 5 Myr and 9 Myr to 7 mm flux density 
would be negligible. The average SFR over this region 
in about the past 3 Myr (before O stars enter the W- 
R phase), based on inferred total mass of young star 
population from 7 mm flux density, is thus estimated as 
4 MQ/yr (Salpeter IMF of 1-100 Mq), in agreement with 
the radio recombination line study by jRodrigucz-Rico et 



of radio SMR ( 


Huang et al. 1994 |Muxlow et al. 


1994 


Kronberg et al. 


2000 l-'enech et al. 


mm. 



5.2. Molecular Gas and Star Formation Efficiency 

Our 7 mm continuum images give a picture of where 
the young star-forming regions are in M82, particularly 
the compact 7 mm sources, which are likely very young 



embedded clusters. What is the spatial relation of these 
sources to the molecular gas? Is there evidence in the 
spatial distribution for how the star formation was trig- 
gered? 

We can estimate M^tan the mass of new stars formed, 
from Ljfi, Mgas, the total mass of gas, from icoi 
and the mass of dense [Mdense] > lO"* cm^^) 

gas from Lhcn- The dense gas fraction is therefore 
Mdense gas/Mgas and the Star formation efficiency, the 
fraction of total gas converted to newly formed stars, is 
Mstar/{Mgtar + Mgas)- This is the Star formation effi- 
ciency for the current generation of stars (i.e., on time 
scales of a few Myr). We convert our observed Sjmm 
in lower resolution map to Lqb in the 2.7" beam, take 
that Lqb to be i/ij, and find Mstar from Mstar = 
1.5 X 10-3 LiR Mq which assumes an O star domi- 
nated population with Salpeter IMF between 1-100 Mq. 
Since we are using the total 7 mm flux density for this 
computation rather than the thermal 7 mm flux den- 
sity, the Lqb and Lm values could be overestimated 
by 30% due to the nonthermal contribution. The full 
translation is uncertain by a factor of 3, mainly corre- 
sponding to the difference in lower mass cutoff in IMF. 
We find the dense gas mass fraction from Ihcn/Ico- 
For the molecular gas mass, we adopt the conversion of 
Mdense gas{H2) = 10 Lhcn Mq {K km s~'^ pc^)-\ and 
M„„ JH2) = 4.8 Lc n Mq {K km s'^ pc^)-^ fromlCao 



& Solomon ( 2004b I . These conversions for molec ular gas 
are unceHain by a factor of 2-3, discussed in |Gao fc:| 
Solomon] poMb) . 

In I'igure [s] we plot the observables Lin/Lco vs. 
Lhcn / Leo ■ The corresponding SFE against the dense 
gas fraction (DGF) are shown on the top and right. The 
data points of M82 represent the results on a 2.7" scale 
(or ^ 46 pc). We show local normal spirals, luminous in- 
frared galaxies (LIRGs), ultra-luminous infrar ed galaxies 



Gao et al. 



ULIR Gs) , and high-z starburst galaxies from 
2007| spirals, luminous infrared galaxies (LIRGs), ultra- 
luminous infrared galaxies ( ULIRGs), and high-z star- 



burst galaxies from Gao et al., (2007) on the same plot. 
The SFE and DGF track m IVlSTroughly as they do in 
all sources, although the size scales of the observations 
can be very different. All M82 data area lie under the 
SFE:DGF=1:1 line, which suggests that less than 50% 
of the molecular gas has been turned into stars. Most 
of the data points lie above the SFE:DGF = 1:10 line, 
implying a limit of ~ 10% SFE for dense molecular gas. 

6. REMARKS ON COMPACT NON-THERMAL SOURCES 

We have detected 5 non-thermal sources with fiux den- 
sities higher than 5 a. Their fiux densities at wavelengths 
between 7 mm and 18 cm in 0.2" x 0.2" beams are listed 
in Table 2\ and plotted in Figure [3] We overlay on the 
measurccTSED the best fit to a model of a single power 
law with free- free absorption. 

The model results are shown in Table E] About one 
third of the compact 7 mm sources we detected (Ta- 
ble 1) are non-thermal. These sources have been known 
for some time and are believed to be SNR (Kronb erg fc:| 
Sramek 19851. However, they do not decline with time 
as expected, so it has been sug gested that they may 
be w ind-driven bubbles instead (Seaquist & Stankovic 
2007|. 
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6.1. 43.992+59.70: Potential AGN m M82? 

The brightest 7 mm continuum compact source is 
43.992+59.70, which has been suggested to be the nu- 
cleus of M82 (Kronberg fc Wilkinson^l975j ) . On the basis 
of European VL_t5i iNetwork and VLA images, Wills et aT] 
( |1999a| suggest this source could be a candidate for the 
AGiN in M82 based on the jet-like feature, although they 
detected only the northern part. The study of two radio 
OH satellite lines also suggests the OH features of this 
source could be associated with a circumnuclear disk or 
torus dSeaquist et"ar]|1997[) . 4 3.992-H59.70 is also within 
3" of the 2 ^m peak ot M82 ( [Lester et ar||1990[ ) that is 
presumably in the cent er of the mass d istribution. In ad- 
dition, [Ne H] line (Achtermann fc Lacy 1995), and OH 



maser kinematics (Seaquist et al 
investigation (|Seaquist et af 



ni i97), molecular gas 
1997 ~and radio recombi- 



ico et al.|2004h also sug- 
the dynamical center of 



nation line results (Rodriguez-Kii 
gest that this source is close to the dyr 
M82, which supports the AGN scenario. 

The suggestion that 43.992+59.70 is an AGN can be 
also tested by loo king for variabili t y. Lo ng-term radio 



flux monitoring by Kronberg et al. ( 2000| suggests that 
the source has a stable flux density, with variability < 
5%, which is the uncertainty in the absolute flux scale 
calibration. Our 7 mm flux density is about 63% of 
the archival centimeter-wave SED model's prediction. It 
is possible that 43.992+59.70 is variable at the current 
epoch, or that the spectrum is not power law to 7 mm. 

A discrepancy from the AGN picture is that Chandra 
high resolution X-ray images do not show 43.992+59.70 
as a strong source. Also, another argument against the 
AGN nature of 43.992+59.70 can be based on the recent 
high resolution VLBI and MERL IN images of this source 
( [McDonald et al.||2001| [Fenech et al. 2008), where it has 
a resolved round structure, consistent with an SNR as 
suggested by its spectral index (a 
However 



al. 19941. 



—0.5) (Huang et 
the measured 7 mm flux density 
appears to miss a signiflcant amount of flux density from 
the SED model of SNR, and it is hard to match the 7 mm 
flux density into the SNR picture. 

This bright non-thermal source also coincid es with a 
compact Pa a and 1.644 /xm [Fell] emitter. iSeaq uist] 
& Stankovic ( 2007 1 modeled the radio emission of non- 
thermal sources m M 82, including this one, as wind- 
driven bubbles. This scenario hypothesizes that non- 
thermal synchrotron emitting particles are produced at 
the shock between the wi nd of young star clu sters and 



hot bubble gas. However, Fenech et al. (20081 measure 



the expension velocity of this source to be at 2700 km/s 
which is ab out 50 times higher than th e expect wind- 
driven shell ( Seaquist fc Stankovic|[2007 l, and it is more 
closer to the tree-expending yiNR. 

In summary, these three models involving AGN, SNR, 
and non-thermal bubble of a young cluster could all ex- 
plain the strong 7 mm emission of 43.992+59.70. De- 
spite some suggestive dynamical evidences, the morphol- 
ogy and expansion revealed by the VLBI and MERLIN 
work strongly indicate that the source is an SNR. How- 
ever, the SNR scenario could not explain the decline of 
flux density at 7 mm from the expectation value of SED 
modeling. The nature of this source might be still open 
to debate. 



6.2. Intermediate Black Hole Candidate M82 X-1 
(CXO M82 J095550.2+ 694047) 



iMatsumoto et al. (2001) reported 9 X-ray sources in 
82. 'i'he source CXITM82 J0 95550.2+694047. or M8 2 
X-1 (09''55'"50.2^ +69°40'47"; [Matsumoto et al.[|200l| , 
has been categorized as an ultraluminous X-ray source 
(ULX) and a 54 mHz quasi-period ic oscillatio n has been 
found in X-ray by Strohmayer & Mushotzky' (2003 1 and 
a possible 62 day X-ray luminosity variation by Kaaret et| 



aT] ([20061). It is suggested to host an intermediate mass 



Elack hole ( Ebisuzaki et al.[|2001 Portegies Zwart et al 



counterpart ot M82 X-1 (Matsumoto et al.|2001 


Kording 


et al. 


2005 


Kaaret et al. 


2006|. 'i'he source 41 


30+59.6 



part. We did not detect a signiflcant radio counterpart of 
M82 X-1 at 7 mm to a Umit of 0.7 mJy/bm. The SEDs of 
the nearby 7 mm sources 40.933+58.98 and 41.628+58.01 
are clearly thermal and are also not candidates for the 
radio counterpart of M82 X-1. 

6.3. Variable Compact Sources 
The flux density of nonthermal sour ce 41.942+57.62 



(41.9 +58 in Kronberg et al. [l985 and Kronberg et al 
2000 1 at 7 mm is signiflcantly below the po wer law model 



Srame 



31 



This source is known to vary (Kronberg fc 
^ 8.4 - 8.8% 



1985); its flux density declines at 
at '--^ 5 GHz, as has been shown and discussed by 
several authors dKronb erg fc Sramek| |1992| l Allen,il999| 
Kronberg et al. 2000; Muxlow et arp 005; Fe nech et al] 
20081). The cm data listed in Table IJ and plotted in Fig- 
were obtained before 1996 plllen 199 9]. Based 
e epoch-corrected 1993 values of Allen fc Kron- 




berg[ (1998 ) for the flux density of 41.942+57.62, and 
tor a spectral index of -0.72 (from best-flt result without 
7 mm data point), we estimate that the flux density of 
41.942+57.62 would have been 8 mJy at 7 mm in 1993, 
three times higher than the observed value. This is con- 
si stent with a de cline of 8.4%/year found at 5 GHz. 

flrst 



Singer et al. 



reported the supernova 
SN 2004am "(UT3S:4U3TT +69:40:48. 1), which was classi - 
fied spectroscopically as Type II (Mattila et al. |2004 l. 
Beswick et al. (2004) report 3-a upper limits tor this 
source ot 0.34 inJy/bm (8.4 GHz) and 1 mJy/bm (15 
GHz) on 2003 November 14. The null detection in their 
MERLIN observations at 5 GHz on March 9 2004 gave 
an upper limit of 0.18 mJy/bm (3 a), and suggested that 
the SN2004am should reach the radio emission peak with 
a significant delay (up to about 1 year) . We can place a 
3 fj upper limit of STmm = 0.6 mJy/bm to the 7 mm flux 
density of SN2004am in our observation of 2005 April 03. 

There are at least 7 known sources for which the 
7 mm flux density predicted from their cm-wave flux 
densities should be larger than 1 mJy (5 — la). 
They are 40.68+551, 43.82+628, 45.41+637, 45.48+647, 
46.18+677, 46.52+639, and 47.37+680. The last four of 
these are at > 20" from the phase center where the noise 
level is a factor of 2 higher than elsewhere and these may 
be down in the noise. Sources 40.68+551 and 43.82+628 
are, however, within the central primary beam in a 
low-noise region. Both 40.68+551 and 43.82+628 are 
SNR. Source 40.68+551 has minor variability in the long 
term radio continuum flux experiment of ,Kronberg et 
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aL] ( |2000[ ), but no variability of 43.82+628 has been re- 
ported, and this may be the first indication. 

For three of the 7 mm sources (39.659+55.73, 
43.992+59.70, and 45.153+61.34) the 7 mm flux den- 
sities are stronger at the epoch of our observation (2005) 
with < 2" beam than at year 2000 -2001 in the obser- 
vation of Rodriguez-Rico et al. (2004) with a 0.6" beam. 
The flux density ditterences are not conclusive 2 a 
of the C-array), but suggest that these sources may be 
varying. The latter two sources are believed to be SNR, 
based on their radio SEDs. Neither of them vary strongly 
at 6 c r n in t he 12-year monitoring program of Kr onberg| 
eTall pOOOl between 1981-1992. Source 39.659+55.73 
is known to vary at 6 cm, fluctuating 20% continuously 
during this time period. But we have characterized this 
source as an H ii region based on its rising spectrum, and 
it therefore should not vary. This source might be asso- 
ciated with a wind-dri ven bubble, as suggested for the 
non-thermal sources bylSeaquist & Stankovic (20071. 



7. CONCLUSIONS 

We present 7 mm VLA continuum maps of the star- 
burst in M82 at two resolutions, 0.2" (B-configuration) 
and 2" (D-configuration) . We recover essentially all of 
the 7 mm flux density of M82 in the lower resolution 
D-array map. An estimated 2/3 of this flux density is 
thermal free-free emission from the starburst. We detect 
less than 10% of the total 7 mm flux density in the high 
resolution B-array map, which is sensitive to only bright, 
compact emission sources. 

We identify 14 compact sources greater than 
0.5 mjy/bm {5a) in the high resolution image, and have 
constructed SEDs using archival VLA data. Nine of the 
compact sources are compact Hii regions, with emission 
measures EM > 10^ cm^^pc, sizes of < 3 pc, and rms 
densities < nl >^^^> 3000. These compact sources are 
located in regions of heavy visual obscuration, with no 
associated near-IR continuum counterparts. This could 
explain the lack of W-R features in M82. Our detection 
limit is Niyc ^ 0.6 x 10^^ s~^, so the smallest of these 
clusters has ~ 100 O stars, and the largest, ~ 1200 O 
stars {Niyc ~ lO^'^ sec^^). The detected compact ther- 
mal sources account for ~ 6% of the total 7 mm free-free 
flux density. Since the maps detect only the most lumi- 
nous compact Hii regions, this ratio should be higher if 
fainter ones are taken into account. We therefore infer 



that the compact Hii region phase in M82 must last at 
least 6% of the Hli region lifetime, '--^l Myr or longer. 

In addition to the compact Hli regions, the high res- 
olution map reveals extended (2", or 30 pc) regions of 
free-free emission with EM > 10^ cm~^pc. The com- 
parison of intensity (EM) and flux (Ni^yc) suggests that 
these regions are likely to be non-spherical in geometry. 
We suggest two possibilities, either these regions indicate 
regions of distributed star formation, with many small 
H II regions possibly reflecting enhanced star formation 
at the ends of the bar in M82, or these regions are sheet- 
like structures, possibly caused by winds or supernovae. 

We have compared lower resolution 7 mm VLA maps 
to maps of molecular gas tracers and dust. The good cor- 
relation between the 7 mm emission and Spitzer IRAC 
Channel 4 map (dominated by 8 fiia PAH features) sug- 

gests a relation NLyc - 2.2 x lO^i s-^{^J\^), for 
intensities /s/im > 0.06 mJy/pc^. This indicates that 
the PAH feature is a good tracer of star formation in 
M82, although we obtain a different relation between 
N^yc and the 8 fim emission than previous work. 

The correlation of 7 mm emission and molecular gas 
tracers (CO and HCN) at 2" scales is excellent, indicating 
that the relation of star formation and HCN em ission 
found on global scales by ( Gao Sz Solomon||2004a I holds 
down to CMC scales. We estimate that the tib'H m M82 
on CMC scales is ~ 1-10%, similar to Galactic values, 
and that the regions of highest SFE also have the highest 
fraction of dense (n^/^ ?J 10'* cm~^) gas. 
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Fig. 1. — M82 radio continuum maps at 7 mm. (a) Naturally-weighted map from the VLA B-configuration data. The color scheme is 
indicated by the color wedge on the top, with flux density values in mjy/beam. The synthesised beam is 0.19" X 0.15", PA= —10°. This 
map is sensitive to structures up to ~ 3.5". The peak flux is 5.2 mjy/bm. (b) Robustly-weighted map made from the D-configuration 

data. Synthesized beam is 1.99" X 1.47", PA = 17° beam. The 6max is ~ 20". Contours are half-integral powers of 2 X 1.12 mJy/bm. 
Dashed box shows the region displayed in (a). 
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Fig. 2. — SED of 9 thermal sources found in 7 mm high resolution (0.2") map. Line represents the best fit single power law SED except 
in the cases of 40.933+58.98 and 42.193+59.13 where the line indicates free-free self-absorption. 
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Fig. 3. — SED of 5 non-thermal sources found in 7 mm continuum high resolution (0.2") map. Line represents power-law synchrotron 
radiation with free-free absorption. 




Fig. 4. — Compact radio sources overlaid on false color VRI image from Svitz er SINGS project. The magenta contours at integers X 
0.5 mjy/beam show 2 cm emission reproduced from VLA archive data (AM671 jMcDonald et al.||2002| at resolution ~ 0.2". The green 
contours at integers X 0.8 mJy/beam present the 7 mm sources shown in FigurelTC! 
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TABLE 1 

Compact 7 mm Continuum Sources: Positions and Fluxes 



Source^ 


Other Name'' 


a(J2000)= 


<5(J2000)<= 
69°40'+ {") 


cPeafc :d 
'-'7mm ' 


Sjmm 


Size° 




(B1950) 


(A99/M02) 


09''55'"+ {") 


(mJy beam"^) 


(mJy) 


II II O 




39.277+54.32 


39.29+54.2 


48.051 


40.64 


0.9 


1.4+0.4 


0.31, 0.25, 


60 


39.659+55.73 


39.68+55.6 


48.427 


41.97 


3.6 


4.6+0.3 


0.21, 0.18, 


170 


/in riQQ 1 c^o no 
4U.yoo+5o.9o 


40. 94+08. o 


/in cif\A 
49.094 


45.20 


0.6 


2.3+0.6 


C\ A'~7 n Q 


70 


41.145+56.29 


41.17+56.2 


49.904 


42.51 


0.8 


3.1+0.6 


0.38, 0.30, 


170 


41.628+58.01 


41.64+58.4 


50.375 


44.22 


1.2 


3.5+0.5 


0.33, 0.26, 





41.942+57.62 


41.95+57.5 


50.686 


43.77 


2.6 


2.1+0.2 


0.18, 0.13, 


180 


42.072+58.52 


42.58+58.4 


50.817 


44.69 


1.2 


2.2+0.3 


0.28, 0.19, 


10 


42.193+59.13 


42.21+59.2 


50.945 


45.32 


1.4 


2.6+0.3 


0.26, 0.21, 


130 


42.549+58.13 


42.56+58.0 


51.291 


44.29 


0.9 


9.2+1.3 


0.68, 0.44, 


140 


42.673+58.32 


42.69+58.2 


51.414 


44.47 


0.6 


9.3+1.5 


0.82, 0.52, 


140 


43.173+58.43 


43.18+58.3 


51.908 


44.56 


0.7 


1.2+0.3 


0.25, 0.19, 


90 


43.295+59.27 


43.31+59.2 


52.028 


45.42 


1.6 


2.0+0.3 


0.22, 0.17, 


170 


43.992+59.70 


44.01+59.6 


52.723 


45.78 


5.2 


5.4+0.2 


0.20, 0.16, 


180 


45.153+61.34 


45.17+61.2 


53.875 


47.41 


1.1 


1.9+0.3 


0.23, 0.21, 


170 



^ Conv ention of Kro nbcig et al. (1985). 

^ A99: [Allenl ( [l999| ; M02: [McDonald et al"] | |2002| |. 

^ The typical uncertainty of position within the 7 mm map is ~ 10 - 20 mas for strong sources, and 30 - 40 mas for 
weaker sources (flux density ^ 1 mjy/bm). The absolute astrometry accuracy is estimated to be ~ O.l", limited by 
the atmospheric phase stability during the observation. 

^ ^7 mm S7 mm' pcak flux dcusity and total integrated flux density of 7 mm continuum, respectively, with beam 
as 0.19 X 0.15 , PA— -10°. The rms uncertainty of each measurement is 0.13 mJy/bm. The systematic error of flux 
density is expected to be < 5%. 

^ Convolved size. Single Gaussian is assumed for each source; using AlPS task IMFIT. The typical uncertainty of source 
size (convolved) is < 10% except source 42.69+58.2 and 42.56+58.0 in which the uncertainty is ~ 20%. 
^ Peaks in emission complex. 



TABLE 2 
Radio SEDs of Compact Sources 



Source^ 






'5'3.6cm 


S2cm 


•S*! .3cm 


Sjmm 


Reference'^ 


(B1950) 


(mJy) 


(mJy) 




(mJy) 




(mJy) 




(mJy) 




(mJy) 




39.277+54.32 




1.7 ± 


0.2 






1.5 ± 


0.2 






1.5 ± 0.5 


M02 


39.659+55.73 




1.5 ± 


0.2 






3.2 ± 


0.2 






8.0 ± 0.6 


M02 


40.933+58.98 




1.9 ± 


0.2 


2.9 ± 


0.1 


2.8 ± 


0.3 


2.2 ± 


0.3 


2.5 ± 0.4 


A99 


41.145+56.29 




2.0 ± 


0.3 






2.8 ± 


0.2 






5.1 ± 0.6 


M02 


41.628+58.01 




1.0 ± 


0.2 






1.8 ± 


0.2 






6.1 ± 0.5 


M02 


41.942+57.62'' 


79.0 ± 4.0 


40.0 ± 


0.9 


26.8 ± 


0.5 


17.9 ± 


1.7 


12.0 ± 


0.6 


2.4 ± 0.3 


AK98 


42.072+58.52 




1.2 ± 


0.4 






1.7 ± 


0.2 






3.6 ± 0.6 


M02 


42.193+59.13 


2.5 ± 1.0 


4.1 ± 


0.5 


4.3 ± 


0.3 


4.5 ± 


0.3 


4.1 ± 


0.3 


4.0 ± 0.5 


AK98 


42.549+58.13 




2.7 ± 


0.4 






1.8 ± 


0.2 






1.7 ± 0.2== 


M02 


42.673+58.32 




1.0 ± 


0.4 






1.8 ± 


0.2 






1.1 ± 0.2<= 


M02 


43.173+58.43 


8.5 ± 1.0 


5.2 ± 


0.2 


3.8 ± 


0.2 


2.4 ± 


0.3 


1.8 ± 


0.3 


1.6 ± 0.5 


AK98 


43.295+59.27 


19.5 ± 1.0 


10.3 ± 


0.2 


8.0 ± 


0.1 


5.5 ± 


0.3 


4.1 ± 


0.4 


2.5 ± 0.5 


AK98 


43.992+59.70 


9.5 ± 1.0 


21.8 ± 


0.9 


20.2 ± 


0.1 


15.8 ± 


0.3 


12.5 ± 


0.3 


5.4 ± 0.4 


AK98 


45.153+61.34 


12.5 ± 1.0 


7.9 ± 


0.3 


5.7 ± 


0.1 


3.9 ± 


0.3 


2.5 ± 


0.3 


2.2 ± 0.5 


AK98 



^ Convention of Kronberg et al. (1985) 

^ S\: Total integrated flux density of radio continuum at wavelength A with 0.2" to 0.3" beam. See ^for details. 
^ The flux densities at 18 — 1.3 cm are adopted from AK98: Allen 8z Kronberg (1998; and references therein); A99: 
Allen et al (1999); M02: McDonald et al. (2002). 



The source is a known strong variable. The flux density densities measured at diffe rent epochs have been corrected 
to the epoch of 1993.4 using a decline rate of 8.8% (see [Allen Sz Kronberg||l998[ for details), except 7 mm flux 
density, which is measured in 2005 in this work. 

^ Source is in emission complex. The peak flux density at 7 mm obtained with IMFIT is used. This could be an 
underestimate if the source is resolved. 
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TABLE 3 
SED Models of Compact Sources 



Source 


Type^ 


Model*^ 


a" 




EM^ 


NLyc" 


No 


Lob 


(B1950) 








(GHz) 


(10"^ pc cm,-^) 


10^" 




10» Lq 


39.277+54.32 


H 11 


PL 


-0.1 


< 5 


< 90 


22 


220 


0.4 


39.659+55.73 


H II 


PL 


0.8 


> 43 


> 8000 


116 


1160 


2.3 


Art OQQ 1 t^O OQ 


XJ TT 
Xl II 






4.9 


85 


36 


360 


0.7 


41.145+56.29 


Hii 


PL 


0.5 


> 43 


> 8000 


74 


740 


1.5 


41.628+58.01 


Hii 


PL 


1.0 


> 43 


> 8000 


88 


880 


1.8 


42.072+58.52 


Hii 


PL 


0.6 


> 43 


> 8000 


52 


520 


1.0 


42.193+59.13 


Hii 


FFSA 




2.5 


20 


58 


580 


1.2 


42.549+58.13 


Hii 


PL 


-0.2 


< 5 


< 90 


25 


250 


0.5 


42.673+58.32 


Hii 


PL 


-0.1 


< 5 


< 90 


16 


160 


0.3 


41.942+57.62' 


SNR 


PL+FFA 


-0.7/-0.9' 












43.173+58.43 


SNR 


PL+FFA 


-0.7 


0.9 


2.7 








43.295+59.27 


SNR 


PL+FFA 


-0.6 












43.992+59.70 


SNR 


PL+FFA 


-0.6 


2.2 


6.6 








45.153+61.34 


SNR 


PL+FFA 


-0.7 


1.0 


3.0 









^ Based on the previous cm results from [Alien & KronbergHl998( ; |Allen| p999[ ; [McDonald et al.H2002[ . 
^ PL: single power law; FFSA: free-free emission with self-absorption; FFA: foreground free-free absorption. 

a: spectral index, S,^ oc i/" . i^t-: turnover frequency. 
^ EM OC 1. ^ 10,000 K is assumed. 

Assuming optically thin frec-frcc emission. For optically thick sources, the numbers would be the lower limits. 
^ Variable source, a — —0.7: with epoch corrected (1993.4) data at wavelengths at 18 - 1.3 cm. a — —0.9: with 
all radio continuum data, including flux density at 7 mm of epoch 2005. See |6.3| 
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